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Abstract 

We show that the Schrodinger equation in phase space proposed by 
Torres-Vega and Frederick is canonical in the sense that it is a natural con¬ 
sequence of the extended Weyl calculus obtained by letting the Heisenberg 
group act on functions (or half-densities) defined on phase space. This al¬ 
lows us, in passing, to solve rigorously the TF equation for all quadratic 
Hamiltonians. 


1 Introduction 


The theory of phase-space Schrodinger equation 
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ih—^f(x,p,t) = H(x + ihdp, —ihd x )'^>(x,p,t) (1) 
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proposed by Torres-Vega and Frederick mmm (TF) has deservedly received 
much attention these last years (see for instance .01? E! ID. Eli USD- TF 
arrive at their equation by using a generalized version of the Husimi transform 
in which coherent states are used as pseudo-transition matrix elements via a 
certain kernel K cs . (We take the opportunity to point out that this kernel 
has been extensively studied in the mathematical literature under the names of 
“wavepacket transform” or “FBI transform”; see e.g. II3I)- We intend to show 
that TF’s choice corresponds to a natural extension of Weyl calculus, and is 
therefore in a sense canonical. 
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Notations 

We will use the collective notations x = (aq,..., x n ), p = (pi, ...,p n ), and z = 
(x,p). We denote by a the canonical symplectic form on the phase space R^ ra = 
1" x K": 

cr(z, z') = (z') T Jz , J = 

if z = (x,p), z' = {x'p'). We denote by Sp{n) the real symplectic group; it 
consists of all linear automorphisms S of R? n such that S T JS = SJS T = J. 
<S(R n ) is the Schwartz space of rapidly decreasing functions on R” and its dual 
S'(M. n ) the space of tempered distributions. The scalar product of two n -vectors 
x and p is denoted by xp, and we will write d z = (d x ,d p ) where 



2 Heisenberg—Weyl Operators on Phase Space 

Let us recall the definition of the Heisenberg-Weyl operators from the “Schrodingerian” 
point of view (see e.g. |S|). The flow (/ t ) determined by the Hamiltonian 
H Zo {z) = cr(z,zo) is given by ft(z) = z + tzo hence the time-one mapping fi 
is just the phase-space translation T(zo) ■ z i—> z + zo; the action of that 
translation on functions T = T(^) is defined by “push-forward”; 

T(zoMz) = #(z-zo). 

The quantized version of the Hamiltonian H Zo is the operator H Zo = a(z, zq) 
with z = (x, — ihd x ), and the solution of the corresponding Schrodinger equation 

Q, 

ih— = a(z,z 0 )ip , ip(x,0)=ip 0 (x) 


0 I 
-I 0 


is the function 


ip{x, t) = exp 




ip 0 {x-tx 0 ). 


The value of if) a point x at time t = 1 is denoted by T(zo)iJj(x): 

f(zo)il>(x) = x - a 0 .); (3) 

T(zq) is by definition the Heisenberg-Weyl operator associated to the translation 
T(z 0 ) (see |5j). So much for the traditional point of view. Let us shortly review 
the notion of phase as studied in our paper [5] . Consider a Lagrangian manifold 
Yq in R? 71 ; we assume for simplicity that Vq is simply connected (this restriction, 
which can be alleviated by passing to the universal covering, is of no consequence 
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for the argument that follows). Choosing, once for all, a base point z in Vg we 
define the phase of Vg by 



where the integral is calculated along an arbitrary path leading from z to 2 . 
This defines <p 0 as a function Vg —» K such that dip 0 (z) = pdx if z = (x,p) (Vg 
being Lagrangian, the integral only depends on the homotopy class with fixed 
endpoints of 7 (z, z ) and since Vg is in addition simply connected there is just 
one homotopy class). Let now H be some arbitrary Hamiltonian function; its 
flow (ft) takes Vg to a new Lagrangian manifold V™ after time t. The phase of 
that manifold is (with obvious notations) 

<p(z,t) = <p 0 (f- t (z)) + J pdx-Hdt 

where the integral is calculated along the arc of phase space trajectory leading 
from the point f-t(z) £ Vg at time 0 to the point z £ V" at time t. Choose 
now for H = H zo ; a straightforward calculation (see 0 , Prop. 10) yields 

t 2 

<fi(z, t ) = tp 0 (z - z 0 ) + tp 0 x - —p 0 x 0 - 

We next consider, as is customary in geometric quantization (see Email), a 
“waveform” on Vg, i.e. an expression of the type 

* 0 ( 3 ) = Jp(z) 

where mo(z) is an integer related to the Maslov index and p a function > 0 
onVg (more accurately, one should view v / p as a “half-density” on Vg but we 
can ignore this point here because half-densities and functions transform the 
same way under translations). The action of a Hamiltonian flow (f t ) on U/g is 
given by ft'f'o(z) = \H( z,t ) where 

^(z,t)=e^ t \ m ^ t \f t )^(z) 

the difference m(z, t)—mo(z) is a count of the number of caustic points appearing 
when the flow deforms Vg into V" = /t(Vg). If H = H Zo the manifold Vg is 
translated in phase space without any “bending” so that there will appear no 
new caustics hence m(z,t ) = mo(z) and we thus have 

*(z,t) = f(tzo)*o(z) (4) 

where the operator T(zo) is now defined by 

f(zo)^o(z) = e^ PoX -i poX ^T(z 0 )^ 0 (z). (5) 
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A straightforward calculation shows that the partial differential equation 
satisfied by ^(z,t) is 

d'H 

ih—(z,t) = E V(z,t) , ty{t = 0) = (6) 

at 

where E is the operator E = — pox + iKzod z ;the latter can be rewritten in terms 
of the symplectic form er as 

E = a{x 0 ,po;x + ihdp, -ihd x ) 

so that m is just TF’s phase-space Schrodinger equation m for the Hamiltonian 

//„>• 

The whole argument can actually be reversed: using for instance the methods 
of characteristics one checks that the solution of equation m is precisely 0 
so that we conclude that the imposition of the TF quantization rules m is 
equivalent to the extension 0 of Heisenberg-Weyl operators to phase-space. 

3 Weyl Calculus on Phase Space 

The process called “Weyl quantization” associates to a suitable function (or 
“symbol”) a = a(z) the operator a w = A defined by 

M{x) = ( 2 hiT f a{z 0 )f(z 0 )i>{x)d 2n zo ( 7 ) 

where a is the “alternative Weyl symbol”, that is the symplectic Fourier trans¬ 
form T a a of a: 

d(z) = T a a(z) = (^) n | e-l^a(z')d 2n z'. 

The discussion of Section |3 suggests that we might now be able to make A 
to act, not only on function of x, but also on functions of 0 by defining 

M*) = ( 2 hi)" J ~a(z 0 )T(z 0 )*(z)d 2n z 0 (8) 

where T now is a function of z = (x,p) and T(zq)^'(z) is given by ©; this 
expression makes perfectly sense provided of course, that a and fk belong to 
some adequate spaces; we will assume this is the case. It turns out that this 
reinterpretation of Weyl operators again leads to the TF quantization rules Gi¬ 
lt is sufficient to prove this in the case of one degree of freedom; dropping indices 
let A' and P be the Weyl operators with symbols x and p; a simple calculation 
of (symplectic) Fourier transforms yields 

X*{z) = ihj6{x 0 )®6'(po)T(zo)*(z)d 2n z 0 

= -ih J 6(zo)(f*(z) - %{z))d 2n z 0 
= {x + ih-§P)^(z) 
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and we thus have X = x + by a similar calculation we get P = 

and we have thus recovered the TF rules 0 . 

4 Exact Solutions for Quadratic Hamiltonians 

The extended Weyl calculus constructed in previous Section allows us to solve 
exactly the TF equation when the Hamiltonian H is a homogeneous quadratic 
polynomial in the variables Xj,pk, in which case the Hamiltonian flow consists 
of symplectic matrices St- 

One of the main features (and probably one of the most attractive!) of 
the usual Weyl correspondence a <—*■ a w = A is the property of metaplectic 
covariance: if we replace abyao S ' -1 where S is any element of Sp(n) then : 

(a o S -1 )™ = Sa w S~ 1 (9) 

where S is any of the two metaplectic operators ±S associated to S (see 0, 
Section 6.3; for the properties of Mp(n) see 0, Chapter 6). The proof of 
property 0 is based on the formula 

Sf(z 0 )S~ 1 = T(Szq) (10) 

where one uses the traditional definition 0 definition of T(zq), hence if itlol) still 
holds when we redefine T(zq) by 0, then metaplectic covariance © will also 
hold for our enlarged Weyl calculus. Now, Mehlig and Wilkinson have shown 
m that if det(S — I) ^ 0 then 

S = ±k J f(Sz 0 )f(-z 0 )d 2n z 0 (11) 

where k is a constant given by 

k=(^ K ) n i v V\det(S-I)\ (12) 

and v is a kind of “Maslov index”. We have shown [B] that every S £ Mp(n) 
can be written as the product of two operators of the type o above (which 
are only defined for det(S l — I) ^ 0); we may thus extend every S G Mp(n) to 
an operator defined on functions of z = {x,p)\ the natural domain of S is then 
5(R^ ra ), and can be extended to L 2 (M. 2n ) and to 5 , (M? n ) by continuity. Having 
done this it is immediate to check, using formulae m that m, and hence also 
©, remain true in our extended Weyl calculus. 

Let us now return to the flow (St); denote by (St) the unique one-parameter 
group in Mp(n ) determined by (St); viewing the St as extended Weyl operators 
we automatically have 

ih—St = H(x + ihdp , — ihd x )St 

hence *1 '(z, t) = SfT o(z) is the solution of the TF equation 0 for the considered 
Hamiltonian. 
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5 Discussion and Concluding Remarks 

We have established that the Torres-Vega and Frederick’s equation fits as part of 
an extended Weyl calculus, which was itself motivated by considerations from 
geometric quantization. There remains to find the physical interpretation of 
that equation. A clue might be the fact that it is possible to rewrite (a variant 
of) the TF equation in terms of the star-product, familiar from deformation 
quantization (this fact seems to have been noticed in this and the fact that 
the star-product is used to express the time-evolution of the Wigner function 
via a quantum-Liouville equation might be an indication that the TF equation 
plays in quantum mechanics a role similar to that of Hamilton’s equation in 
classical mechanics. 

We will come back to this important question in a near future. 
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